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ABSTRACT 

In the solar photosphere, many properties of coronal holes (CHs) are not 
known, especially vector magnetic fields. Using observations from Hinode, we 
investigate vector magnetic fields, current densities and current helicities in two 
CHs and compare them with two normal quiet regions (QRs) for the first time. 
We find that, in the CHs and QRs, the areas where large current helicities are 
located are mainly co-spatial with strong vertical and horizontal field elements 
both in shape and location. In the CHs, horizontal magnetic fields, inclination 
angles, current densities and current helicities are larger than those in the QRs. 
The mean vertical current density and current helicity, averaged over all the 
observed areas including the CHs and QRs, are approximately 0.008 A m~^ 
and 0.005 m~^, respectively. The mean current density in magnetic fiux 
concentrations where the vertical fields are stronger than 100 G is as large as 
0.012 ± 0.001 A m~^, consistent with that in the fiare productive active regions. 
Our results imply that the magnetic fields, especially the strong fields, both in 
the CHs and QRs are nonpotential. 

Subject headings: Sun: activity — Sun: magnetic fields — Sun: photosphere 
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1. 



Introduction 



Coronal holes (CHs) are low density and temperature regions in the solar atmosphere 
(Munro & Withbroe 1972). They appear as dark areas observed with X-ray or EUV lines. 
CHs are always classified into three categories according to their locations and lifetimes: 
polar, nonpolar (isolated), and transient (Harvey & Recely 2002). In previous studies, many 
properties of CHs and their relationship with magnetic fields have been investigated by 
many authors, such as temperature variation (Wilhelm 2006; Zhang et al. 2007), element 
abundance (Laming & Feldman 2003), magnetic field evolution (Yang et al. 2009a, b), and 
magnetic field structures (Meunier 2005; Zhang et al. 2006; Tian et al. 2008). 

In CHs, magnetic fields are predominated by one polarity and open magnetic lines are 
concentrated (Bohhn 1977). Plasma escapes along the open magnetic flux, giving rise to 
fast solar wind (Krieger et al. 1973; Tu et al. 2005). CHs differ from the normal quiet 
Sun regions mainly through the difference of magnetic structures, i.e., the field lines are 
mainly closed in the quiet Sun and open above CHs (Altschuler et al. 1972). Hot gas is 
able to escape along open field lines but is trapped in closed loops. The normal quiet Sun 
appears brighter due to the radiation of the trapped gas. However, the magnetic fields are 
not exclusively unipolar in CHs and consequently CHs should contain also locally closed 
coronal loops besides the open flux (Levine 1977). Wiegelmann & Solanki (2004) computed 
some properties of coronal loops. They found that high and long closed loops are extremely 
rare, whereas short and low-lying loops are almost as abundant in CHs as in the quiet Sun. 
This result suggests an explanation for the relatively strong chromospheric and transition 
region emission (many low- lying, short loops), but the weak coronal emission (few high and 
long loops) in CHs. 

Vector magnetic fields are very significant since they can provide us plentiful information, 
such as electric current and current helicity. In a common consensus, the free energy is 
stored in the stressed, nonpotential magnetic fields. Electric current and current helicity are 
two parameters which are always used to characterize magnetic nonpotentiality (Moreton 
& Severny 1968; Abramenko et al. 1996; Wang 1996; Zhao et al. 2009; Su et al. 2009). 
Electric current density J is calculated as 



where /iq is the magnetic permeability in vacuum (47rxl0 ^ G m A By definition, current 
helicity density he is derived as 



The larger the deviation from the potential fields is, the more the energies that can power 




VxS, 



(1) 



K = BVxB. 



(2) 
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solar activity are. Therefore, electric current and current helicity also serve as measures of 
the productivity of solar activity. 

Electric current and current helicity have been intensively studied in active regions 
(ARs), while quantitative study is quite rare in CHs. This is mainly due to the relative 
weak fields in CHs and the former instruments are not capable of achieving reliable vector 
field observations. Fortunately, Hinode (Kosugi et al. 2007) provides us wonderful space- 
based measurements which are vastly superior to previous data in resolution and sensitivity. 
The direct purpose of this study is to compare vector magnetic fields, current densities and 
current helicities in CHs and quiet regions (QRs) with Hinode observations. We describe the 
observations and data reduction in Section 2. The results are presented in Section 3, and 
the conclusions and discussion in Section 4. 



2. Observations and data reduction 

The data used here were carried out with the Spectro-Polarimeter (SP; Lites et al. 2001) 
in the Solar Optical Telescope (SOT; Ichimoto et al. 2008; Shimizu et al. 2008; Suematsu et 
al. 2008; Tsuneta et al. 2008) instrument aboard Hinode. The SP provides Stokes /, Q, U, 
and V profiles of two Fe lines at 630.15 nm and 630.25 nm with wavelength sampling of 21.6 
mA in four modes (normal, fast, dynamics and deep maps). We adopt the observations of 
two CHs (named CHI and CH2 here) and two non-coronal hole QRs (QRl and QR2) taken 
in the fast map mode. For each raster slit, the integrated exposure time is 3.2 s and the pixel 
sampling along the slit is 0".32. The scan is in the east-west direction with a scanning step 
of 0".30. CHI was observed from 11:42:25 UT to 12:45:36 UT on November 30, 2008, and 
QRl from 17:34:06 UT to 18:37:18 UT on April 12, 2007. Both of them have a field-of-view 
(FOV) of 300" X 162". They are located almost at the disk center. CH2 was taken from 
01:36:29 UT to 02:08:49 UT on November 12, 2007, and QR2 from 11:35:09 UT to 12:21:47 
UT on April 13, 2007. Both of them have a FOV of 151"x 162". They are centered at about 
10° from the disk center. 

The SP data are calibrated and inverted at the Community Spectro-polarimetric Analy- 
sis Center (CSAC; [http : //www . csac . hao . ucar . e du/ ) . Vector magnetic fields are retrieved 
using the inversion techniques based on the assumption of the Milne-Eddington atmosphere 
model (Kubo et al. 2007; Yokoyama, T. 2010, in preparation). In the inversion, a non-linear 
least-squares fitting technique is used to fit analytical Stokes profiles to the observed profiles. 
Values of 36 parameters are returned from the inversion, including the three components of 
magnetic field (field strength B, field inclination 7, field azimuth 0), the stray light fraction 
a, and so on. In the vector field measurements based on the Zeeman Effect, there exists a 
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180° ambiguity in determining the field azimuth. Various algorithms have been developed 
to resolve the ambiguity, but it is also difficult to make a complete removal (Metcalf et al. 
2006; Semel & Skumanich 1998). As pointed out by Metcalf et al. (2006), "the methods 
which minimize some measure of the vertical current density in conjunction with minimiz- 
ing an approximation for the fields' divergence show the most promise." In this study, we 
use the improved Nonpotential Magnetic Field Calculation (NPFC) method developed by 
Georgoulis (2005) to disambiguate the azimuth angles. 

The vector magnetic field B is shown by longitudinal field strength Bcos^ and transverse 
field strength Bsinj without involving filling factor, where B is the intrinsic strength of the 
magnetic field and 7 the inclination angle with respect to the linc-of-sight direction. The 
transverse field can be decomposed into B^—Bsin'ycos(j) and Bj^=Bsin'ysin(j) in the image 
plane, where is the field azimuth angle. Then we transform the vector magnetic field at 
each pixel to heliographic components {B^, By, B^) according to the formulae given by Gary 
& Hagyard (1990). Besides, geometric mapping of the magnetograms in the image plane into 
the heliographic plane is also performed. Vector magnetic fields in the photosphere allow us 
to compute the vertical component of electric current J in the photospheric plane where 
z = 0. Jzi i.e., 

is calculated by finite differences. Then we obtain the current helicity density he as, 

hc{x, y,z^O)^ fJ'oBzJz\z=o. (4) 

Here, he is just part of the current helicity density since only the vertical current density can 
be calculated. 

We determine the noise level of vertical fields by calculating la standard deviation of 
weak fields within intranetwork regions and find it is about 4 G. The noise level of horizontal 
fields determined in the similar way is about 35 G. For current density, the noise level (0.005 
A m~^) is estimated from the standard deviation of them in the pixels where the horizontal 
fields are weaker than 35 G, similar to the method of de La Beaujardiere et al. (1993), Wang 
et al. (1996) and Wheatland (2000). For current hehcity density, the noise level (0.004 G^ 
m~^) is calculated from the standard deviation of them over the pixels where either the 
vertical or the horizontal fields are weaker than their noise levels. 



3. Results 



CHI and QRl are located at the disk center, while CH2 and QR2 have nearly the 
same heliocentric angle. We investigate the vector magnetic fields, current densities and 
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current helicities of each pair of targets, respectively, and the unsigned quantities are used 
throughout. 



3.1. Comparison between CHI and QRl 

CHI is dominated by the positive polarity as shown in the upper panel in Figure 1. The 
dash-dotted curve encloses CHI area the pixels within which are considered only. We plot 
the probability density functions (PDFs) of the vertical fields in CHI and QRl in Figure 
2 (upper panel). Not unexpectedly, CHI fields have a generally imbalanced PDF between 
the positive and negative polarities. The positive fields, which are the predominant polarity, 
have a higher PDF than the negative. We define the magnetic flux imbalance parameter p 
as 

3) ... 

where j) is the vertical magnetic flux at pixel (i, j). As presented in the lower panel, the 
flux imbalance in CHI is about 0.65, much higher than that in QRl (0.25). Figures 3 and 4 
display the appearance of horizontal magnetic fields and the distribution of derived vertical 
current helicities in CHI and QRl, respectively. When we only consider the areas where 
unsigned vertical magnetic fields are stronger than 4 G and horizontal fields are stronger 
than 35 G, the mean vertical magnetic field, horizontal field, current density and current 
helicity in CHI are 33.27 ± 0.01 G, 96.64 ± 0.07 G, 0.00888 ± 0.00001 A m-^, and 0.00492 ± 
0.00001 G^ m-\ while those in QRl are 39.44 ± 0.01 G, 89.20 ± 0.06 G, 0.00839 ± 0.00001 
A m~^, and 0.00558 ± 0.00001 G^ m~^, respectively. Here, the uncertainties of the means 
are calculated according to Equation (4.14) in Bevington & Robinson (2003). 

We examine the distribution of current hehcities (Figure 4) in CHI and QRl and find 
that the large current helicity concentrations are mainly co-spatial with the strong vertical 
and horizontal fleld elements (Figures 1 and 3). The region highlighted with square "1" in 
Figure 4 is investigated in detail. At the center of this region, large current hehcities are 
distributed. In Figures 1—3, the contours are vertical magnetic fields at ± 100 G levels. At 
the places where large current helicities are located, the vertical and horizontal magnetic 
fields (squares "1" in Figures 1 and 3) arc much stronger than the surrounding areas. We 
take 100 G as the general separation of strong and weak magnetic fields. The mean vertical 
field strength within the contour is 584.1 ± 0.2 G and the mean horizontal field strength 
is 162.3 ± 1.3 G. The corresponding average unsigned current density and current helicity 
are 0.0130 ± 0.0002 A m'^ and 0.0962 ± 0.0002 G^ m-^ respectively. In QRl, we select 
another sub-region (square "2" in Figure 4) where large current helicities are concentrated. 
Similar to region "1" in CHI, the areas with large current helicities also correspond with 
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strong vertical and horizontal fields. The mean vertical field, horizontal field, electric current 
and current helicity within the contour curve in region "2" are 593.8 ± 0.1 G, 176.7 ±1.3 
G, 0.0131 ± 0.0002 A m-^, and 0.0988 ± 0.0001 G^ m"!. 

To examine whether these magnetic properties change with vertical magnetic field 
strength, we calculate the mean values of horizontal magnetic fields, field inclination angles, 
current densities, and current helicities above some thresholds of unsigned vertical magnetic 
fields in CHI and QRl (Figure 5). With the increase of threshold from weak vertical fields 
to about 100 G, the mean value of horizontal fields in CHI increases sharply from 97 G to 
145 G, as shown in Figure 5a. When the vertical fields exceed 100 G, the horizontal fields 
exhibit a slight increase to the level of 180 G. The relationship between vertical fields and 
horizontal fields in QRl is similar to that in CHI except for about 15 G weakness. The mean 
inchnation angle of the magnetic lines in CHI is about 75° for all the fields (Figure 5b). Then 
it decreases quickly to about 38° with the increase of vertical fields in the range from weak 
to about 100 G. When the vertical fields are stronger than 100 G. the inclination angles 
gradually decrease to less than 8°. While the mean angles in QRl are about 5° smaller that 
those in CHI. In CHI, the current densities increase with the vertical fields and then stay 
at the level of 0.014 A m~^, as presented in Figure 5c. While the average current density in 
QRl is about 0.012 A m~^. The difference between them is relatively significant considering 
their uncertainties. The current helicities in CHI reveal a general trend of increase when the 
vertical fields increase from weak fields to kilogauss (Figure 5d). We notice that the current 
helicities in QRl exhibit similar variation trends to those in CHI with smaller means. 



3.2. Comparison between CH2 and QR2 

CH2 is an enhanced network field region with predominance of negative polarity (top 
left panel in Figure 6) and QR2 a normal quiet Sun with mixed-polarity (top right). We 
also examine the PDFs of their vertical magnetic fields. Similar to those in CHI and QRl, 
the PDF in CH2 is clearly asymmetric while the PDF in QR2 shows balanced positive and 
negative magnetic fields (upper panel in Figure 7). Compared with CHI in Figure 2, CH2 
shows a much more serious imbalance between the opposite polarities (lower panel in Figure 
7). The flux imbalance p in CH2 is as high as —0.8 and that in QR2 is only about —0.1. The 
mean vertical magnetic field, horizontal field, current density and current helicity in CH2 
are 41.21 ± 0.01 G, 86.94 ± 0.07 G, 0.00808 ± 0.00001 A m^^ and 0.00543 ± 0.00001 G^ 
m-\ while those in QR2 are 33.34 ± 0.01 G, 84.42 ± 0.08 G, 0.00796 ± 0.00001 A m^^ and 
0.00429 ± 0.00001 G^ m-\ respectively. 

In Figure 6, the bottom panels are corresponding current helicity distributions in CH2 
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and QR2 scaled between ± 0.03 A m~^. Prom their general appearance, we can see that 
large current helicities are also mainly co-spatial with strong magnetic fields both in shape 
and location. We select two areas (highlighted with windows "3" and "4" ) containing large 
current helicity patches. When we take the pixels with vertical fields stronger than 100 G 
into consideration, the mean vertical field strength in region "3" is as strong as 733.1 ± 
0.2 G, and the mean horizontal field, current density and current helicity are 152.8 ± 1.6 
G, 0.0118 ± 0.0002 A m-^, and 0.1153 ± 0.0002 G^ m'^, respectively. While these four 
parameters in region "4" are 511.3 ± 0.3 G, 109.0 ± 2.2 G, 0.0121 ± 0.0003 A m-^, and 
0.0824 ± 0.0003 G^ m-\ respectively. 

We also compare the variations of horizontal fields, inclination angles, current densities 
and current hehcities versus the thresholds of unsigned vertical magnetic fields in CH2 and 
QR2. The horizontal fields increase both in CH2 and QR2 before the vertical field threshold 
reaches 100 G (Figure 8a). Then the horizontal fields in CH2 increase slowly from 120 G 
to 150 G, while they mainly maintain at the 115 G level in QR2. The inclination angles in 
both CH2 and QR2 decrease with the increase of vertical fields in a similar trend (Figure 
8b). The inclination angles finally approach to 5°. The trend of current density variation 
in CH2 is consistent with that in QR2 (Figure 8c). They are similar to those in CHI and 
QRl shown in Figure 5c. When the threshold of vertical fields is several hundred Gauss, the 
mean of current densities in CH2 is larger than that in QR2. The current helicities in both 
CH2 and QR2 reveal a general trend of increase when the vertical fields increase from weak 
fields to kilo-Gauss though the rise slows down (Figure 8d). Generally, the current densities 
in CH2 are larger than those in QR2. 



3.3. Comparison between the CHs and QRs 

Figure 5 indicates the means of the horizontal magnetic fields, inclination angles, current 
densities, and current hehcities in CHI are larger than those in QRl. While Figure 8 exhibits 
that in CH2, only the horizontal magnetic fields and current hehcities are stronger than QR2, 
and the inchnation angles and current densities in them are the same for the stronger fields. 
To enlarge sample, we combine two CHs (CHI and CH2) and two QRs (QRl and QR2), 
respectively. We compare the parameters between the CHs and QRs (Figure 9), and find 
that horizontal magnetic fields, inclination angles, current densities and current helicities in 
the CHs are larger than those in the QRs. 



-8- 



4. Conclusions and discussion 

Using the Hinode/SP data, we investigate vector magnetic fields, current densities and 
current helicities in two CHs (CHI and CH2), and compare them with two QRs (QRl and 
QR2). To our knowledge, this comparison has not been done using vector field measurements 
before. We find that in the areas both in the CHs and QRs where large current helicities 
are concentrated, there are strong vertical and horizontal field elements and they are mainly 
co-spatial with each other in shape and location. In the CHs, horizontal magnetic fields, 
inclination angles, current densities and current helicities are larger than those in the QRs. 
The mean values and their uncertainties of five parameters in the areas where the unsigned 
vertical magnetic fields are stronger than 4 G and horizontal fields than 35 G are listed in 
Table 1. Averaged over the observed CHs and QRs, the means of vertical magnetic fields, 
horizontal fields, inclination angles, current densities and current helicities are approximately 
37 G, 90 G, 75°, 0.008 A m-^ and 0.005 G^ respectively. 

The mean vertical magnetic strength in CH2 is 41 G, stronger than that in CHI (33 
G). The reason is that CH2 is located at an enhanced network region. According to the 
MDI Synoptic Chart, the location of CH2 is ~250° (Carrington longitude). At the same 
site, there was a small mature AR on 17 October 2007, which has been studied by Murray 
et al. (2010). We think that the dispersion of the AR's negative polarity leads to the strong 
fields in CH2. The average vertical magnetic field in QR2 is 33 G, somewhat weaker than 
that in QRl (39 G). It may be caused by the magnetic field fiuctuation of the quiet Sun 
in the scale of FOV of our data. As revealed in our results, the CHs are dominated by one 
polarity, while magnetic fiuxes in the QRs are generally balanced in sign. 

There is also another way to estimate whether the difference or similarity is significant, 
i.e. using the uncertainties of the vector field parameters from the inversion. But this method 
is not available, because "some of the errors of the inversion parameters are wrong and the 
problem can not be fixed any time soon," as pointed out by Lites (private communication). 
"However, the values for the inversion parameters themselves are OK." 

Flaring ARs carry more current than simple ones which are closer to potential. Ac- 
cording to Mcintosh (1990), the majority of regions are simple. As presented in our results, 
both in the CHs and QRs, large current helicity areas are almost located in the regions 
with strong vertical and horizontal magnetic fields. The mean current density in magnetic 
concentrations where the vertical fields are stronger than 100 G is as large as 0.012 ± 0.001 
A m~^, consistent with that (~ 0.01 A m~^) in the flare productive ARs (de La Beaujardiere 
et al. 1993; Leka et al. 1993; Wang et al. 1994; Wang et al. 1996; Zhang 2001; Deng 
et al. 2001; Leka & Barnes 2003; Liu et al. 2007; Li et al. 2009). These results imply 
that the photospheric magnetic flelds, especially the strong flelds, in the CHs and QRs are 
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nonpotential. 

Since only two pairs of CHs and QRs are investigated here, more observations will be 
taken into consideration to check the similarities and differences between CHs and QRs. 
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Table 1. Mean values and their uncertainties of five parameters in the areas with unsigned 
vertical magnetic fields stronger than 4 G and horizontal fields stronger than 35 G. 





Vertical field 


Horizontal field 


Inclination 


Current density 


Current helicity 




(G) 


(G) 


(Degree) 


(A m-2) 


(G2 m-i) 


CHI 


33.27 ± 0.01 


96.64 ± 0.07 


77.14 ± 0.03 


0.00888 ± 0.00001 


0.00492 ± 0.00001 


QRl 


39.44 ± 0.01 


89.20 ± 0.06 


74.98 ± 0.02 


0.00839 ± 0.00001 


0.00558 ± 0.00001 


CH2 


41.21 ± 0.01 


86.94 ± 0.07 


73.41 ± 0.03 


0.00808 ± 0.00001 


0.00543 ± 0.00001 


QR2 


33.34 ± 0.01 


84.42 ± 0.08 


74.93 ± 0.03 


0.00796 ± 0.00001 


0.00429 ± 0.00001 
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Coronal Hole 30-Nov-2008 1 1 :42:25— 12:45:36 UT 




Quiet Sun 12-Apr-2007 17:34:06— 18:37:18 UT 




Solar X (arosec) 



Fig. 1. — Overview of the vertical magnetic fields in CHI (upper panel) and QRl (lower 
panel). The dash-dotted curve delineates the CH boundary. Squares "1" and "2" outline 
two sub-regions that are investigated in detail. The contours in two squares are vertical 
magnetic fields at ± 100 G levels. 
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Fig. 2. — Probability density function of the vertical magnetic fields (upper panel) and flux 
imbalance (lower panel) in CHI and QRl. 
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Fig. 3. — Appearance of the horizontal magnetic fields corresponding to Figure 1. The 
dash-dotted curve, squares, and contours have the same meanings as those in Figure 1. 
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Coronal Hole 30-Nov-2008 1 1 :42:25.3— 12:45:36.4 UT 




Solar X (arosec) 



Fig. 4. — Distribution of the vertical current helicities in CHI and QRl. The dash-dotted 
curve, squares, and contours have the same meanings as those in Figure 1. 
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Fig. 5. — Variations of the horizontal magnetic fields (a), field inclinations (b), current 
densities (c), and current helicities (d) versus the thresholds of unsigned vertical magnetic 
fields in CHI and QRl. Each error bar represents the standard deviation of the mean value. 
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Fig. 6. — Vertical magnetograms in CH2 (top left) and QR2 (top right) and corresponding 
current helicity distributions (bottom panels). Windows "3" and "4" highlight two areas 
that are studied in detail. 
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Fig. 7. — Probability density function of the vertical magnetic fields (upper panel) and fiux 
imbalance (lower panel) in CH2 and QR2. 
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Fig. 8.— Similar to Figure 5 but for CH2 and QR2. 
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Fig. 9.— Similar to Figure 5 but for the CHs (CHI and CH2) and QRs (QRl and QR2). 



